Quantum Andreev map: A paradigm of quantum chaos in superconductivity 
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We introduce quantum maps with particle-hole conversion (Andreev reflection) and particle-hole 
symmetry, which exhibit the same excitation gap as quantum dots in the proximity to a super- 
conductor. Computationally, the Andreev maps are much more efficient than billiard models of 
quantum dots. This makes it possible to test analytical predictions of random-matrix theory and 
semiclassical chaos that were previously out of reach of computer simulations. We have observed 
the universal distribution of the excitation gap for large Lyapunov exponent and the logarithmic 
reduction of the gap when the Ehrenfest time becomes comparable to the quasiparticle dwell time. 
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One-dimensional (ID) quantum mechanical models 
with a chaotic classical limit were first studied by Casati, 
Chirikov, Ford, and Izrailev in 1979 1]. These mod- 
els have since developed into one of two paradigms of 
quantum chaos 0, U The other paradigm is the 2D 
billiard of irregular shape 0. Two is the lowest num- 
ber of dimensions for non-integrable (chaotic) dynamics 
in autonomous systems, since a single constant of mo- 
tion is sufficient for integrability in ID. The ID models 
get around this constraint through a periodically time- 
dependent external force ("kick"), which eliminates the 
energy as constant of motion — but still conserves the 
quasi-energy (analogously to quasi-momentum conserva- 
tion in a periodic lattice). The two paradigms share a 
common set of phenomena in the fields of quantum chaos 
and localization [HE 00. 

The combination of chaos and superconductivity pro- 
duces an entirely new phenomenology, notably the ap- 
pearance of an excitation gap as a signature of quan- 
tum chaos :9]- The paradigm common to most of the 
literature is the 2D billiard connected to a superconduc- 
tor ^Jj, introduced under the name "Andreev billiard" 
in Ref. [ll|. The name refers to the Andreev reflection 
which occurs at the interface with the superconductor, 
where an electron at energy e above the Fermi level is 
converted into a hole at energy e below it. 

From the point of view of computational efficiency, 
compact quantum maps such as the kicked rotator [lj 
(a particle confined to a circle and driven periodically 
in time with a strength that depends on its position) 
are much more powerful than 2D models such as bil- 
liards. Indeed, there exists a highly efficient diagonaliza- 
tion technique that works only for maps [l2T | . The lack of 
a ID map for quantum chaos with superconductivity has 
hindered the numerical test of a variety of analytical pre- 
dictions E3, EM EH E3, EM E3| Most notably, nu- 
merical efforts have not been able to distinguish the con- 
flicting predictions of random-matrix theory (RMT) 
and the semiclassical Bohr-Sommcrfcld (BS) quantiza- 
tion: RMT predicts an excitation gap at the Thouless 
energy while BS gives an exponentially vanishing den- 
sity of states without a true gap. Recent analytical work 



13 0) El Ell has predicted that diffraction creates 
a gap in the BS density of states at the inverse Ehrenfest 
time. This has never been seen in computer simulations, 
because the Ehrenfest time scales logarithmically with 
the system size and is usually far too small to play a 
role. For these reasons there is a real need for something 
like a "quantum Andreev map". Does it exist? And 
if it does, can it be simulated more efficiently than the 
Andreev billiard? These are the issues addressed in this 
paper. 

We show how to construct quantum Andreev maps 
out of any conventional quantum map (not necessarily 
chaotic), in much the same way as any normal billiard 
can be turned into an Andreev billiard by coupling it to 
a superconductor. The construction is guided by the clas- 
sical electron and hole dynamics on the Poincare surface 
of section of an Andreev billiard. The Andreev kicked 
rotator is a particular example of such an Andreev map. 
We certify that it possesses the phenomenology of the 
Andreev billiards and search for predictions of RMT and 
semiclassics. We leave for future investigations the appli- 
cation of the Andreev map to other kicked models (pos- 
sibly with a different phenomenology), such as the kicked 
top 0] and the Fermi-Ulam model |21| . 

A quantum map is represented by the Floquet op- 
erator F, which gives the stroboscopic time evolution 
u(ptq) — F p u(0) of an initial wave function u(0). (We set 
the stroboscopic period tq = 1 in most equations.) The 
unitary M x M matrix F has eigenvalues exp(— ie m ), 
with the quasi-energies e m £ (— 7r, tt) (measured in units 
of Ti/tq). This describes particle excitations in a nor- 
mal metal. We also need hole excitations. A particle 
excitation with energy e m (measured relatively to the 
Fermi level) is identical to a hole excitation with en- 
ergy — E m . This means that hole excitations in a nor- 
mal metal have Floquet operator F* and wave function 
v (p) = (F*)Pv(0). 

Particles and holes are coupled by connecting the nor- 
mal metal via a lead to a superconducting reservoir. On 
a Poincare surface of section, the lead is represented 
by a spatially localized region in which Andreev reflec- 
tion converts electrons into holes and vice versa, with 
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phase shift —i. (A weak energy dependence of this phase 
shift is ignored for simplicity, but can be accounted for 
straightforwardly.) Analogously, for the quantum An- 
dreev map we assume that Andreev reflections occur 
whenever an excitation ends up in a certain subspace 
of Hilbert space. This subspace Tlx, riz, . . . njv consists of 
N out of M states in a chosen basis and corresponds to 
a lead with N propagating channels. The N x M matrix 
P projects onto the lead. Its elements are P nm = 1 if 
?7i = n G {n 1 ,7i 2 , ■ ■ .tt-at} and P nm = otherwise. The 
dwell time of a quasiparticle excitation in the normal 
metal is Tdwell — M/N, equal to the mean time between 
Andreev reflections. The fact that Andreev reflections 
only occur at multiples of the stroboscopic time to is 
technically convenient, and should be physically irrele- 
vant for Tdwell > TQ. 

Putting all this together we construct the quantum 
Andreev map from the matrix product 

T-v( F ° V-( 1 - PTP - IPTP ) (1) 
* ~ 1 { F* ) ' ' - { -iP T P 1 - P T P ) ■ (1 > 

(The superscript "T" indicates the transpose of a ma- 
trix.) The particle-hole wave function W = (u,v) evolves 
in time as &(p) — ^^(O). The Floquet operator can be 
symmetrized (without changing its eigenvalues) by the 
unitary transformation T — > p~ 1 / 2 J 7 'P 1 / 2 1 with 

,pi/ 2 _ / 1 - (1 - W2)P T P ^V2P T P \ 
1 ~\ -i±V2P T P 1 - (1 - ±V2)P T P J ■ 

(2) 

In order to establish the correspondence of the ID 
quantum Andreev maps to 2D Andreev billiards, we ex- 
amine the spectral properties of the map. The Floquet 
operator T possesses a particle-hole symmetry which en- 
tails that its 2M eigenvalues exp(— ie m ) come in inverse 
pairs. This symmetry is the analogue of the particle-hole 
symmetry in Andreev billiards, in which excitation en- 
ergies ±e occur symmetrically around the Fermi level. 
To avoid double-counting of levels, we restrict the quasi- 
energy to the interval (0,7r). The excitation spectrum 
of particles and holes consists of the M quasi-energies 
in this interval, and the mean level spacing tt/M is 
twice as small as the level spacing 6 — 2ir/M for par- 
ticles and holes separately. The energy scale for the 
proximity-induced excitation gap is the Thouless energy 
E T = NS/An = N/2M = l/(2T dwcll ). 

The quantization condition dct(jF — e~ l£ ) = can be 
written equivalently as 

det[l + 5(e)5*(-e)] = 0, (3) 

in terms of the N x TV scattering matrix [22| 

S(s) = P[e~ l£ - F(l~ P T P)}- 1 FP T . (4) 

Eq. © for the Andreev map has the same form as for 
the Andreev billiard |23| , but there S is given in terms of 
a Hamiltonian Hq instead of a Floquet operator F. The 



two approaches become entirely equivalent in the context 
of RMT, when Ho is chosen from one of the Gaussian en- 
sembles and F is chosen from one of the circular ensem- 
bles [24[- They are also equivalent in the semiclassical 
limit, when the billiard can be represented by a Poincare 
map which can be quantized approximately |25j| . 

In the mean-field limit M > N > 1, RMT predicts a 
hard gap in the excitation spectrum of size -Ermt = 7i?T 
(with 7 = 2- 3 / 2 ( v / 5- If/' 2 = 0.60), and above the gap a 
square-root singularity in the density of states p(e) = 
Tr-iA" 3 / 2 ^ - £rm T ) 1/2 (with A = QMS^I^S) 0. 
Sample-to-samplc fluctuations of the lowest excitation 
energy £n around the mean-field gap have been calculated 
in Refs. |13lll4l |. A universal probability distribution was 
predicted for the rescaled energy x = (eq — £'r.mt)/A. 
While the mean-field prediction of RMT has been tested 
numerically in an Andreev billiard , the numerical er- 
ror bars are too large to extract the predicted universal 
gap fluctuations. 

To demonstrate the efficiency of the quantum Andreev 
maps, we specialize to the quantum kicked rotator. The 
Floquet operator is 

{Mo d 2 \ ( KI \ 

xexp (4t^) ! (5) 

with To the moment of inertia of the particle and K 
the (dimcnsionless) kicking strength. The particle moves 
freely along the circle for half a period, is then kicked 
with a strength K cos 8, and proceeds freely for another 
half period. The transition from classical to quantum 
behavior is governed by the effective Planck constant 
h c g = Htq I Iq . Since we would like to compare the kicked 
rotator to a chaotic billiard, without localization, we fol- 
low the usual procedure of quantizing phase space on 
the torus 9,p £ (0, 2ir), rather than on a cylinder, with 
p = —ihcffd/dB, the dimensionless angular momentum 
|3 . For h cB = 2ir/M, with integer M, the Floquet oper- 
ator is an M x M unitary symmetric matrix. In angular 
momentum representation it has elements 

F kk ,=e-W 2M X k2+k ' 2 XUQtf) kk ,, (6a) 
U kk , =M- 1 / 2 e ( 2 «/M)fefc' ; (6b) 

Ofefe' =4 fc ' 6 - (lMK/27r)cos(27rfe/M) . (6c) 

Upon increasing K the classical dynamics varies from 
fully integrablc (K = 0) to fully chaotic [K > 7, with 
Lyapunov exponent A w \n(K/2)]. For K < 7 stable and 
unstable motion coexist (a so-called mixed phase space). 

To introduce the Andreev reflection we use a projection 
operator which is diagonal in p-representation, 

, dTd , c / 1 if L < k < L + N- 1, 

( F p )kk>=S*k>x[ otherwise. ^ 

(We checked that similar results are obtained when P is 
diagonal in ^-representation.) The position L of the lead 
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FIG. 1: Main plot: Gap distribution for the Andreev kicked 
rotator with M = 8192, K = 45, and M/N = r dwc n = 10 
(diamonds), 20 (circles), 40 (+), and 50 (x). The solid line 
gives the RMT prediction [l3||. Inset: Density of states for 
the same system. The solid line is the RMT prediction 
The dashed line is a numerical result in the mixed regime 
(M = 8192, K = 1.2, M/N = 10). 
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FIG. 2: Density of states for the Andreev kicked rotator with 
M = 131072, Tdwcii = 5, and K = 14 (solid line), compared 
with the Bohr-Sommerfeld calculation (histogram), and the 
RMT prediction (dashed line). 



to the superconductor is arbitrary. The Floquet opera- 
tor T of the "Andreev kicked rotator" is then obtained 
by inserting Eqs. I© and Q into Eq. We apply 

the symmetrization (|2J) , so T is a unitary symmetric ma- 
trix. The real symmetric matrix + T^) can be diago- 
nalized efficiently with 0(M 2 In M) operations (and not 
0(M 3 ) as with standard methods) by means of the Lanc- 
zos technique, if the multiplication with the matrix U is 
carried out with the help of the Fast-Fourier- Transform 
algorithm [l2j. The eigenvalues cose m uniquely deter- 
mine the quasi-energy e m G (0,7r). 

In the inset to Fig. ^ we show the density of states 
p(e) for system size M — 8192, kicking strength K = 45 
(strongly chaotic dynamics), and several widths N of the 
lead to the superconductor. The density of states has 
been averaged over 250 different positions of the lead. 
The data points fall on top of the RMT prediction Q 
without any adjustable parameter. Reducing the kicking 
strength down to K = 1.2, one enters the regime of mixed 
classical dynamics. We see that the gap disappears, as 
predicted in Ref. (2f|. 

To test RMT beyond the mean-field limit, we study 
the statistical fluctuations of the gap. The main panel of 
Fig. ^ shows the probability distribution of the smallest 
eigenvalue Eq i n the chaotic regime. To improve statis- 
tics we sampled 6000 different positions of the lead. We 
rescaled the energy x = (eq — £rmt)/A. as prescribed by 
Ref. Good agreement is observed with the univer- 

sal scaling distribution [2jJ, again without any adjustable 
parameters. 

It is predicted theoretically that deviations from RMT 
should appear if the Ehrenfest time te = A -1 InAf is no 
longer small compared to Td we ii- For te ^ th wrii the semi- 
classical Bohr-Sommerfeld approximation [!J |2g should 



be valid, with a diffraction induced gap of order Ti/te, 
[Ttij . To search for these deviations from RMT we con- 
sider rotators with smaller kicking strengths (but still in 
the fully chaotic regime), thus smaller Lyapunov expo- 
nent, and much larger M. 

In Fig.|21we show the density of states for M = 131072 
and K = 14. Strong deviations from the RMT predic- 
tion are clearly visible. Also plotted is the result of a BS 
calculation [jg, in which we slightly smoothed the sin- 
gular delta functions. This approximation agrees better 
with the exact result. Most remarkably, it reproduces 
the three distinct peaks in the density of states, which 
now can be identified with trajectories of certain lengths. 
All trajectories with lengths that are odd multiples of 
Tdwcii = 5 contribute to the peak at e/Et — ir, odd mul- 
tiples of 4 contribute at e/Et — 57r/4, and odd multiples 
of 3 contribute at e/Et — 57r/3. 

A systematic reduction of the excitation gap is ob- 
served upon increasing the ratio te /Tdwellj a s shown in 
Fig. |21 The main panel is a semi-logarithmic plot of 
e /E T as a function of M G [2 9 , 2^ for M/N = r dwoU = 
5 and K = 14. Existing theories 0,H(j predict a linear 
initial decrease of £o with InM at fixed Tdweii = M/N. 
We fit the data to the prediction of Vavilov and Larkin 

— ^2— = 1- ° flnM-21n^-A (8) 
-Ermt 2Ardwcii V N J 

We find a — 0.59 and a 1 = 3.95. Once a and a' are ex- 
tracted, no free parameter is left, and the resulting curve, 
shown with a solid line in the inset to Fig. |3 correctly 
reproduces the parametric dependence on Tdwcii = M/N 
at fixed M. As a further check, we tried the slightly 
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FIG. 3: Main plot: Dependence of the mean gap on the 
system size M, for Tdwoii = M/N = 5 and K — 14. Averages 
have been calculated with 400 (for M = 512) to 40 (for 
M > 5 • 10 5 ) different positions of the contacts to the super- 
conductor. The error bars represent the root-mean-square of 
eo. The dashed line is the RMT prediction and the solid line 
is a linear fit to the data points. Inset: Dependence of the 
mean gap on r dwc ii = M/N for K = 14 and M = 524288. 
The dashed line is the RMT prediction and the solid and 
dotted curves are given by Eqs. JHJ and respectively, 
with coefficients extracted from the linear fit in the main plot. 

different expression 

-=^- = l-^—(hxM-a"), (9) 

»RMT ^ATd wc ll 



with a" — a' + 2 In 5. The resulting curve (dotted line 
in the inset to Fig. shows significant deviations from 
the numerical data. We conclude that Eq. (jHJ) gives the 
correct parametric dependence of the Andreev gap. 

A discrepancy remains in the value of the numerical 
coefficients. While the coefficient a' is model dependent, 
the prefactor a is expected to be universal. Our numerics 
gives a = 0.59 ± 0.08, in between the two competing 
predictions a = 0.23 of Ref. 20] and a = 2 of Ref. 

In conclusion, we have constructed a quantum map 
that accounts for the presence of superconductivity. The 
"Andreev kicked rotator" introduced above, has been 
shown to be equivalent to the Andreev billiards stud- 
ied so far. Owing to the fact that it is one-dimensional 
rather than two-dimensional, it is much more efficient 
computationally, which permits to observe two theoreti- 
cal predictions that are currently out of reach of billiard 
simulations: The universal gap fluctuations for large Lya- 
punov exponent and the logarithmic reduction of the gap 
for small Lyapunov exponent. We foresee that the An- 
dreev kicked rotator on a cylinder (instead of on a torus) 
can be an equally effective tool to study the interplay of 
superconductivity and localization. 

We have benefitted from discussions with I. Adagideli 
and J. Tworzydlo. This work was supported by the 
Dutch Science Foundation NWO/FOM. 
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